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Abstract 

We apply twisted boundary conditions to lattice QCD simulations of three-point 
correlation functions in order to access spatial components of hadronic momenta 
different from the integer multiples of 27r/L. We calculate the vector and scalar 
form factors relevant to the K ^ tt semileptonic decay and consider all the possible 
ways of twisting one of the quark lines in the three-point functions. We show that the 
momentum shift produced by the twisted boundary conditions does not introduce 
any additional noise and easily allows to determine within a few percent statistical 
accuracy the form factors at quite small values of the four-momentum transfer, which 
are not accessible when periodic boundary conditions are considered. The use of 
twisted boundary conditions turns out to be crucial for a precise determination of 
the form factor at zero-momentum transfer, when a precise lattice point sufficiently 
close to zero-momentum transfer is not accessible with periodic boundary conditions. 


1 Introduction 


In lattice simulations of QCD the spatial components of the hadronic momenta pj {j = 
1, 2, 3) are quantized. The specific quantized values depend on the choice of the boundary 
conditions (BCs) applied to the quark fields. The most common choice is the use of 
periodic BCs in the spatial directions 


Tp{x + ijL) = Tp{x) , 


( 1 ) 


that leads to 

27r 

Pj=njj^, ( 2 ) 

where the n^-’s are integer numbers. Thus the smallest non-vanishing value of pj is given 
by 27r/L, which depends on the spatial size of the (cubic) lattice {V = L^). For instance a 
current available lattice may have L = 32 a, where a is the lattice spacing, and a~^ ~ 2.5 
GeV leading to 27r/L ~ 0.5 GeV. Such a value may represent a strong limitation of the 
kinematical regions accessible for the investigation of momentum dependent quantities, 
like e.g. form factors. 

In Ref. [1] it was proposed to use twisted BCs for the quark fields 

V^(x -|- CjL) = ip{x) (3) 


which allows to shift the quantized values of pj by an arbitrary amount equal to 9j2Tr/L, 
namely 


Pj 


9 .^ 

L 


+ Pj = (^j-r + 


2tt 

T 


27r 

T 


(4) 


The twisted BCs (3) can be shown [I] to be equivalent to the introduction of a U(l) 
background gauge field coupled to the baryon number and applied to quark fields satis¬ 
fying usual periodic BCs (the Aharonov-Bohm effect). Moreover in Ref. [2] the twisted 
BCs were hrstly implemented in a lattice QCD simulation of two-point correlation func¬ 
tions of pseudoscalar mesons. The energy-momentum dispersion relation was checked 
showing that the momentum shift 2Tr6j jL is a true physical one. 

The aim of the present paper is to explore the use of distinct twisted BCs for different 
fermion species in the calculation of the momentum dependence of lattice three-point 
correlators^. We want both to establish whether twisted BCs are able to provide in 
practice form factors at small values of the momentum transfer not accessible with 
periodic BCs, and to estimate the level of statistical precision that can be achieved. 
The latter is an important point because the introduction of twisted BCs leads to a non- 
negligible increase of the computational time due to the need of producing new inversions 
of the Dirac equation for each quark momentum. We anticipate that the answer to both 

^In this case the equivalent background gauge fields are coupled to the generators in the Cartan 
subalgebra of the flavor group U{N)v commuting with the quark mass matrix. 
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questions is positive; the twisted BCs are able to provide form factors at small values of 
the momentum transfer with a precision comparable to the one attainable with periodic 
BCs. 

In this work we consider the case of the vector and scalar form factors relevant to 
the K ^ TT semileptonic transition, which have been recently investigated using periodic 
BCs by the quenched simulations of Ref. [3], as well as by the n/ = 2 and n/ = 2 + 1 
dynamical flavor simulations of Refs. [4] and [5], respectively. There are two reasons for 
our choice. 

First, the theoretical uncertainty in the determination of the vector form factor at 
zero-momentum transfer, /(O), presently dominates the corresponding uncertainty in the 
extraction of the value of the Cabibbo angle from decays. As shown in Ref. [3], an 
important source of uncertainty for /(O) comes from the precision in the determination 
of the form factor slopes at zero-momentum transfer. Thus an interesting issue is to 
check whether the access to new lattice points at small values of the momentum transfer 
can improve significantly the precision of the determination of /(O). 

Second, the approach of Ref. [3] is characterized by the use of a suitable double ratio 
of three-point correlators, which allows to access in a very precise way the scalar form 
factor fo{q^) at g = 0 corresponding to The precision level of 

the values of /o(<?max) is crucial to achieve the percent accuracy for /(O). However, there 
are cases in which it is not possible to get a lattice point sufficiently close to = 0 using 
periodic BCs. An important example is represented by some of the form factors entering 
hyperon semileptonic decays, like the weak magnetism, the weak electricity, the induced 
scalar and pseudoscalar form factors [6]. Thus, we have carried out an analysis of our 
lattice data for the K ^ n transition, but excluding the very accurate value /o(<?max)- 
We consider this analysis representative of the cases where a precise lattice point close 
to or at zero-momentum transfer is not accessible with periodic BCs. 

We limit ourselves to quenched simulations where the generation of gauge configu¬ 
rations is clearly independent of the BCs applied to the quark fields^. We expect that 
the ability of twisted BCs to provide form factors at small values of the momentum 
transfer do not depend on the use of the quenched approximation, and therefore we are 
confident that our findings hold as well also in case of partially quenched and full QCD 
simulations. 

We obtain the K ^ tt vector and scalar form factors for quite small values of the 
four-momentum transfer, which are not accessible when periodic BCs are considered, 
without introducing any significant additional noise and within a few percent statistical 
accuracy. For completeness we consider all the possible ways of twisting one of the quark 
lines in the three-point functions. 

When the precise lattice point for fo{q^) at q^ = g^ax is not included in the analysis, 
the use of twisted BCs is crucial to allow a determination of /(O) at a few percent 

^It has been recently shown [7, 8] that for many physical processes, including semileptonic decays, 
one can impose twisted BCs on valence quarks and periodic ones for sea quarks, eliminating in this way 
the need for producing new gauge configurations for each quark momentum, since finite-volume effects 
remain exponentially small. 
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statistical level. On the contrary, when the very accurate value /o(9max) i® included in 
the analysis, the impact of twisted BCs on the determination of the K ^ n form factors 
and their slopes at zero momentum transfer turns out to be marginal, while we remind 
it is expensive for the computational time. This result is due to the fact that: i) 

the precision of the lattice points obtained with twisted BCs is not comparable to the 
one that can be achieved through the double ratio method of Ref. [3] at the particular 
kinematical point = q'^ax) which in turn is quite close to = 0 in the simulation of 
Ref. [3]; ii) the use of twisted BCs does not lead to a sufficient improvement of the 
precision with respect to the nearest space-like points obtained with periodic BCs. 

The plan of the paper is as follows. In the next Section we briefly discuss the 
implementation of the twisted BCs for the evaluation of all the propagators required in 
this work. In Section 3 we present the calculation of two- and three-point correlation 
functions having twisted quark lines. In Section 4 we show our results for the scalar 
and vector form factors of the K ^ n transition, while in Section 5 we investigate the 
impact of the twisted BCs on the determination of the slopes of the vector and scalar 
form factors at zero momentum transfer. Finally Section 6 is devoted to our conclusions. 

2 Lattice quark propagators with twisted BCs 

On the lattice, for a given flavor, the quark propagator S{x,0) = (V'(x) ipiO)), where 
(...) indicates the average over gauge field configurations, satisfies the following equation 

^L>(x,y) S’(2/,0) = 4,0 (5) 

y 

where D{x, y) is the Dirac operator whose explicit form depends on the choice of the 
lattice QCD action. In what follows we work with Clover fermions and therefore D(x, y) 
is given explicitly by 

D{x, y) = (amo + 4r) - ^ ^ {^x,y-afi {r - Jy) U^{x) 

y 

+ 4,y+aA + Uliy)'^ - CTyu , (6) 

where Uy{x) is the gauge link, Pyu{x) is the (symmetric) plaquette in the plane 

and we have omitted Dirac and color indices for simplicity. 

When the quark field satisfies the twisted BCs (3), the corresponding quark propa¬ 
gator ^(xjO) = {i^{x) V’(O)) still satisfies Eq. (5) with the same Dirac operator D{x,y) 
but with different BCs. Following Refs. [1, 2] one can redefine the quark field as i’g{x) = 

in order to work always with periodic BCs on the fields. In such a way 
the new quark propagator S^{x,0) = {ipffix) V'e(O)) satisfies the following equation 

J2D^{x,y) S^\y,0) = 5,,o ( 7 ) 

y 
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with a modified Dirac operator D^(x,y) but periodic BCs. The new Dirac operator is 
related to Eq. ( 6 ) by simply rephasing the gauge links 

U,{x) ^ U^{x) ( 8 ) 

with the four-vector 6 given by (0, 6). Note that the plaquette P^y{x) is left invariant by 
the rephasing of the gauge links. In terms of S^{x,y), related to the quark fields 'ijj{x) 
with periodic BCs, the quark propagator S{x, y), corresponding to the quark fields V'(x) 
with twisted BCs, is simply given by 

= . (9) 

3 Two- and three-point correlation functions with twisted 
quark lines 

We are interested in calculating the 7r~ form factors of the weak vector current 

Vfj, = which are defined through the relation 

= f+{Q‘^)iP + P')^^ + f-{Q‘^)iP - p')fM , ( 10 ) 

where = {p — p')‘^. As usual, we express f-{q^) in terms of the so-called scalar form 
factor 

Mq^) = + ( 11 ) 

with /o(0) = /+(0). 

From Eq. (10) the form factors can be expressed as linear combinations of hadronic 
matrix elements of time and spatial components of the weak vector current. The latter 
can be obtained on the lattice by calculating two- and three-point correlation functions 

Cf^nt.,ty,P,p') = ^_^{OAty,y)%{t,,S) ( 12 ) 

S,y 

C^W(t,p-) = E(Oicw(Tx)Ojcw(0))e-*^'-% (13) 

X 

where Oj. = (I 75 U, = dyss are the operators interpolating 7 r“ and mesons, and 

V)j is the renormalized lattice vector current 

= Zv[l + + cvd^ sa^^^u) , (14) 

where Zy is the vector renormalization constant, by and cy are 0 (a)-improvement 
coefficients and the subscript (. refers to the light u (or d) quark. In what follows we 
always use degenerate u and d quarks. 
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Using the completeness relation and taking tx and {ty — tx) large enough, one gets 








{'^ip')\%\K{p)) e 


tx') 


(15) 


C^W(t,p(pO) 


2E 




K (tt) 


(16) 


where Ek = \/+ \p\'^, E^ = ^/M^TWF and = (0|<^i^(7r)(0)l-fi'(7r))- Then 

it follows 


Cf['^itx,ty,p,p') 
C^{tx,P) C^{ty - tx,p') 


tx —*00 

{ty — tx) ^ OO 


{7r{p')\%\K{p)) 

VZrZtt 


(17) 


Consequently the hadronic matrix elements (7r(p')| V)j|-f7(p)) can be obtained from the 
plateaux of the l.h.s. of Eq. (17), once Zr and are separately extracted from the 
large-time behavior of the two-point correlators (16). 

In terms of the strange and light quark propagators Ss{x,y) and S^{x,y) the two- 
point correlator C^^'"\t,p{p')) becomes 

C^W(t,pV)) = (18) 

X 

where 5'^(^)(0, x) = 75S']^^^(x, 0)75. Analogously the three-point correlator C^'^{txity,p,p') 
becomes 


Cf^'^itx^ty^p^p') = ^{Tr[E^{x,Q;ty,p')-i^Ss{Q,x)])e (19) 

X 

where the generalized propagator T,{x,0;ty,p') satisfies the equation [9] 

J2D{x,z) T,{z,0]ty,p') = 755(x,0)e“*^’''^ 5t,^ty ■ (20) 

Z 

Using the twisted propagators, S and S, Eqs. (18)-(20) hold as well taking into ac¬ 
count the corresponding change of the quantized momenta {pj pj). The new two- and 
three-point correlators can be always expressed in terms of quark propagators satisfying 
periodic BCs, namely Eq. (7). For the case of our interest and adopting in what follows 
the convention that the values of the momenta p and p' are always given by multiples 
of 27r/L, we get 

C^\t,‘^{h-ei)+p) = ^(Tr[5|^(x,0)755f(0,x)75])e-*^''" 

X 

C--(t,^(0s-02)+p') = J2('E^l^dH^,0)75S^H0,x)75])e-^E-x ( 2 I) 
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^(03 - Ol) +P, ^(03 - 02) +p’) = 

X 

0; (0, a;)]) 

Cf^\t,,ty, ^(03 - 0-i) +p, ‘fiOs - 02 ) +p') = 

X 

(Tr[Sf'f-^(x, 0; ty,p'hX^ (0, x)]) (22) 

where q = p—p' and the generalized propagator is solution of the modified equation 

Y.dIi{x,z) Sg’j(z,0;t„p') = 754(^>0)e"'^"'" • (23) 

Z 

To help visualizing the notation used for the 0-vectors the three-point correlator relevant 
for the —> TT semileptonic transition is depicted in Fig. 1 . 



Fig. 1: Three-point eorrelation function of the —> vr semileptonic transition with the 

various quark lines twisted by the vectors 6 i, 02 and 63 . 

Note that in the r.h.s. of Eqs. (21)-(23) the exponentials do not contain any 0-vector 
because of Eq. (9). 

4 Results for the K —^ n form factors 

We have generated 100 quenched gauge field configurations on a 24^ x 56 lattice at 
P = 6.20 (corresponding to an inverse lattice spacing equal to a~^ ~ 2.6 GeV), with the 
plaquette gauge action. Using non-perturbatively 0(a)-improved Wilson fermions we 
have chosen quark masses corresponding to one pair of the values adopted in Ref. [3], 
namely k G {0.1339,0.1349}. 

Using K and tt mesons with quark content {kskp} and (k^k^) respectively, two different 
R —> vr correlators have been computed, using both kg < k^ and kg > kg, 
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corresponding to the cases in which the kaon(pion) is heavier than the pion(kaon). Using 
the same combinations of quark masses, also the three-point tt ^ K correlations 
have been calculated. Finally, two non-degenerate K ^ K and two degenerate vr ^ tt 
three-point functions have been evaluated. 

As for the critical hopping parameter, we have adopted the value kc = 0.135820(2) 
found in Ref. [3] using the axial Ward identity, and we have chosen the time insertion of 
the vector current equal to ty = T/2, which allows to average the three-point correlators 
between the left and right halves of the lattice. Using the degenerate K ^ K and vr ^ tt 
transitions, for which the vector form factor at zero-momentum transfer is known to be 
equal to unity because of vector current conservation, the values of the renormalization 
constant Zy and of the 0(a)-improvement parameter appearing in Eq. (14), are 
found to be in good agreement with the corresponding results of Ref. [3]. Finally, we 
adopt for the improvement coefficient cy the same non-perturbative value cy = —0.09 
used in Ref. [3]. Throughout this paper the statistical errors are evaluated using the 
jacknife procedure. 

As shown in Ref. [3] the scalar form factor can be calculated very efficiently 

at = {Mk — M^)^ using a double ratio of three-point correlation functions 

with both mesons at rest, namely 


p . C^nt.,ty,0,0)C^^{t,,ty,0,0) 

■^OyJ'Xi ^y) — 


C^^{t„ty, 0 , 0 ) C^it^ty, 0 , 0 ) 


[foiQ: 


2 

max 


tx) ■ 


2 {Mr + 

4MkM^ ■ ^ ’ 


For / <?max the form factors f+{q‘^) and fo{q^) can be determined from the ma¬ 
trix elements of the time and spatial components of the vector current, (7r|Vo|iF) and 
(7r|Uj|iF) {i = x,y,z), that in turn can be extracted from the plateaux of the l.h.s. of 
Eq. (17) (as discussed in the previous Section). However, such a plain strategy leads to 
a determination of fo{q‘^) with a quite poor precision (see Ref. [3]). Therefore in order 
to achieve a much better accuracy we follow Ref. [3] and introduce a suitable ratio of the 
spatial and time components (7r|Vo|.^) and {7r\Vi\K), normalized by the corresponding 
degenerate K ^ K transition, namely 


^ ^ (ttIUIAT) (ATlUoliF) 

(vrlUoli^) {K\nK) ’ 


(25) 


which allows to access the ratio f-{q^)/f+{q^). 

First of all, as a consistency check of our run, we have evaluated the two- and three- 
point correlation functions (18) and (19) setting all the three 0-vectors to zero and 
adopting the same values of p and p' used in Ref. [3], where a different, larger set of 230 
quenched gauge configurations were employed. The results obtained for both /o(<?max) 
and the scalar and vector form factors at ^ 9max nicely consistent within the 
statistical errors. In what follows we will label these results as “0 = 0”. Then, the two- 
and three-point correlation functions (21) and (22) have been evaluated choosing always 
p = p' = 0 and making three different kinematical choices that correspond to assuming 
non-vanishing only one out of the three 0-vectors of Fig. 1. Defining Q = 27r|0|/L, 

Er = ^-|- Q'^ and we consider 



Kinematics A: 6*1 = 0 / 0 and 6*2 = 6*3 = 0 pK = 2 tt 9/L, p-,^ = 0) 


UiQ^) = 


1 


2 M^ 

/o(g') = f+iq^) 


{7r\Vo\K) - 


1 + 


Q 


{Ek + - {Mk + Ek) 


Ml-M^ Mk + Ek - {Ek - M^)Rio 
= Mk + Ml - 2 Mt,Ek = gmax “ ‘2‘Mt,{Ek - Mk) ; 
Kinematics B: 02 = 0 and 0i = 03 = 0 pK = 6, Pn = 2 tt9/L) 


(26) 


U{q^) = 


1 


2Mk 

/o( 0 ^) = /+( 0 ^) 




1 + 




0 ^ — {Et, + MK)Rio 


Ml - Ml Mk + Ek- (E^ - MK)Rio 
(f = mI + — 2MkEt^ = ~ ‘2‘MK{E-n — M^) ; 


(27) 


Kinematics C: 9^ = 9 ^ 0 and 0i = 02 = 0 (^ = p-^ = 2^9 jV) 

In this case the vector form factor is given by f+{q'^) = {TT\Vi\K)/2Q; however 
a more accurate determination can be obtained by constructing a double ratio 
similar to the one in Eq. (24), but using the spatial components of the weak vector 
current, namely 


Ri{txi ty) — 


Cl^{tx,ty,PK,Pn) Cl^{tx,ty,p^,PK) 

Cl^{tx,ty,PK,PK) Cf"" (tx ,ty,Pn,P7T) 


[Uiqnr (28) 


(ty - tx) ■ 


and 


fo{q^) = f+{q‘^) 


1 + 


1 ( 2Ek 
M l - Ml Ek-eA RiO 


- Ek - Ej, 


q^ = {Ek - Et, 


(29) 


Notice that in kinematics C, by varying |0|, the values of q'^ are always positive in 
the range 0 ^ g^^x- 

We consider two values of |0|, namely |0| = 0.225,0.70. The first value leads to a 
difference + (27r|0|/L)2 — Mk{tt) just exceeding the statistical errors, while the 

second one simply gives rise to a value of which is approximately half of the minimum, 
non-vanishing value attainable for space-like q'^ with periodic BCs. 

For each of the above values we consider two different orientations: 0 = (1,0,0) • |0| 
(asymmetric) and 0 = (1,1,1) • |0|/\/3 (symmetric). In this way we can check whether 
either spatially asymmetric or symmetric momentum shifts lead to different noises. We 
have found no significant difference. 
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In Table 1 we have collected the values of the meson masses, the SU(3)-breaking 
parameter (M^ — M^) and of Q'^ax ™ lattice units that characterize our simulation. We 
remind that for each value of 0 a new inversion of Eq. (7) is required with a computational 
time similar to the one needed for the inversion of Eq. (5) at 0 = 0. 


kg - ke 

aMx 

aiWjr 


22 - 

^ ^max 

0.1339 - 0.1349 

0.3025(20) 

0.2419(24) 

-b0.03299(26) 

0.00367(8) 

0.1349 - 0.1339 

0.3025(20) 

0.3556(15) 

-0.03495(30) 

0.00282(8) 


Table 1: Values of the hopping parameters kg and ki, — M^), and 

The quality of the plateaux used for extracting the vector form factor /+((?^) can 
be appreciated in Eig. 2 for three representative cases corresponding to 0 = 0 and 
6 = (0.7,0,0) in kinematics A and C^. Note that in the latter kinematics he use of the 
double ratio (28) allows to reduce strongly statistical fluctuations. 

In Fig. 3 our results for the form factors f+{q‘^) and obtained at kg = 0.1349 

and ki = 0.1339 for the three kinematics A, B and C as well as for 6 = 0, are reported, 
while in Fig. 4 the relative statistical errors, A/_|_(g^)//_!_ (g^) and A/o((?^)//o((?^), are 
shown. It can be seen that the use of twisted BCs allows to explore the low-g^ region 
without introducing any significant additional noise. 

The statistical error of the results obtained at 0 = 0 quickly decreases as q'^ increases 
in the space-like region, reaching a minimum value of ~ 3% at ~= —0.05 for both 
f+{q^) and fo{q‘^). The results obtained using kinematics A and B are totally consistent 
with each other and the statistical error remains almost constant (~ 4%). Note that 
the precision obtained in kinematics A and B is not better than the one of the nearest 
points obtained at 9 = 0. 

Vice versa a significant improvement in the precision can be achieved using kinematics 
C thanks to the double ratio (28), which allows to determine /+(g^) with a statistical 
error of ~ 0.3%. For the scalar form factor the corresponding accuracy is only ~ 1.5% due 
to the larger fluctuations of the ratio (25), but it is still better than the precision achieved 
with kinematics A and B. The smallest statistical uncertainty (~ 0.07%) remains the 
one at q^ = thanks to the double ratio (24) that involves the time component of 
the weak vector current with all mesons at rest. 

® For sake of simplicity we do not report explicitly the definition of the quantity R+{tx,ty) appearing 
in Fig. 2 for each kinematics. It suffices to say that it is defined in terms of the l.h.s. of Eq. (17) for the 
various components of the weak current which should be combined to determine the vect or form fa ctor 
in the various kinematics [see Eqs. (26-28)]. In kinematics C R+{tx,ty) is defined as Ri{tx, ty) [see 
Eq. (28)]. 


10 





Fig. 2: Time dependence of the quantity R+{tx = t,ty = T/2) providing at large time 
distances the vector form factor (see footnote 3 in the text), calculated with peri¬ 

odic BCs (6 = Oj at afq^ = —0.067 (full squares), and with twisted BCs (6 = (0.7,0,0)J 
at a^q^ = —0.034 in kinematies A (open diamonds) and afq^ = 0.0022 in kinematics 
C (open dots). The plateaux are taken from t/a = 11 to t/a = 17. The values of the 
hopping parameters are kg = 0.1349 and k^ = 0.1339. 


5 Slopes of /_|_(g^) and /o(q^) at = 0 


In this Section we analyze the momentum dependence of the vector and scalar iiT —> vr 
form factors in order to understand the impact of the introduction of the twisted BCs 
in the determination of the form factors as well as of their slopes at zero-momentum 
transfer. As already noted in Ref. [3] the results for f+{q‘^) can be very well described 
by a pole-dominance fit, where the slope agrees well with the inverse of the iF*-meson 
mass square for each combination of the simulated quark masses. The results for fo{q‘^) 
can be parameterized using different functional forms, like the ones adopted in Ref. [3]. 
However, for our purposes it suffices to consider the case of a polar fit for both f+{q‘^) 
and fo{<l'^)- We have checked that our final findings remain unchanged if instead of a 
polar fit a linear or quadratic (?^-dependence for fo{q'^) is considered. 

Thus our results have been parameterized using the following momentum dependen¬ 
cies 


UiQ^) 

Mq'^) 


/(O) 

1 - A+ ’ 

/(O) 

1 - Ao 


(30) 
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2 2 

a q 



Fig. 3: Results for the form factors f+{q^) (a) and fo{q‘^) (h) obtained for various values 
of q^ at kg = 0.1349 and ki = 0.1339. Full squares correspond to 9 = 0, while open 
diamonds, triangles and dots correspond to the three kinematics A, B and C, respectively. 


where /(O) = /+(0) = /o(0), Aq and A+ are fitting parameters. We have determined 
the values of such three parameters through a x^-minimization procedure applied to 
different sets of lattice data, namely: i) data obtained using periodic BCs only (0 = 0); 
ii) addition of the results corresponding to kinematics A and B; iii) data at 0 = 0 
plus those corresponding to kinematics C only; iv) full set of lattice points (0 = 0 plus 
kinematics A, B and C). In Table 2 we have reported the values obtained for /(O), Aq 
and A+ for each choice of the data set. It can be seen that the impact of the lattice 
points corresponding to twisted BCs appears to be marginal and the extracted values 
of /(O), Aq and A+ have almost the same accuracy as the one obtained using simply the 
lattice data calculated with periodic BCs. The usefulness of the twisted BCs is thus 
spoiled by the absence of a significant improvement of the accuracy as well as by the 
non-negligible increase of the computational time required for the inversions of Eq. (7) 
for each choice of 0. 

This finding can be traced back to the following facts: i) the presence of the very 
accurate value of foiq"^) at q‘^ = (/max ™ the data sets considered (such a data point 
is by far the most accurate one [see the rightmost point in Fig. 4(b)] and also the values 
of (Zmax ^re quite small [see Table 1]); and ii) the use of twisted BCs does not lead to 
a sufficient improvement of the precision with respect to the nearest space-like points 
obtained with periodic BCs. 

As mentioned in the Introduction, there are cases of phenomenological interest where 
a lattice point at zero-momentum transfer is not accessible with periodic BCs. An impor- 
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a q 
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a q 
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Fig. 4: As in Fig. 3 hut for the relative statistical errors A/+(( 7 ^)//+(g^) (a) and 

A/o(g2)//o(g2) 


tant example is represented by some of the form factors entering hyperon semileptonic 
decays, like the weak magnetism, the weak electricity, the induced scalar and pseu¬ 
doscalar form factors [ 6 ]. Thus, in order to clarify the role played by the 0 7 ^ 0 lattice 
points, we have repeated the fitting procedure for the same lattice data sets but excluding 
always the very accurate value /o(Q'max)- This analysis is expected to be representative 
of the cases where a precise lattice point close to or at zero-momentum transfer is not 
accessible. Our results are reported in Table 3. 

The following comments are in order: 

• using the data set 6 * = 0 (only periodic BCs) the values of /(O) and of the slopes are 
now more poorly determined. The statistical uncertainties turn out to be ~ 10%, 
~ 40% and ~ 30% for /(O), Aq and A+, respectively. Note that the accuracy 
obtained for /(O) is almost 3 times the precision of the nearest space-like points 

3%); 

• the accuracy improves by a factor of ~ 2 when the data corresponding to kinematics 
A and B are included. Note that the accuracy of /(O) is now comparable to the one 
of the 0 7 ^ 0 points (~ 4%), due to the fact that the latter are quite close to = 0. 
Clearly the ~ 5% precision achieved for /(O) is not enough for phenomenological 
applications to the decay. However such a level of accuracy is highly desirable 
for other observables, like hyperon semileptonic form factors; 

• the inclusion of the quite accurate points obtained in kinematics C (with or without 
those of kinematics A and B) leads to a remarkably good determination of /(O) 
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lattice data set 

/(O) 

Ao/a^ 

A+/a^ 

0 = 0 

0.9907(19) 

3.2(5) 

5.9(8) 

0 = 0 -|- kin. A, B 

0.9920(17) 

2.8(5) 

5.6(7) 

0 = 0 + kin. C 

0.9910(16) 

3.1(5) 

6 . 0 ( 8 ) 

0 = 0 + kin. A, B, C 

0.9921(15) 

2.7(4) 

5.6(7) 


Table 2; Values of the vector form factor at zero-momentum transfer /(O) = /+(0) = 
/o(0) and of the slope parameters Aq and A+ of the polar fits (30), obtained using different 
sets of values o//+(g^) and /o(9^), including always the very accurate value of ffiq^) at 
= 9max- values of the hopping parameters are kg = 0.1349 and = 0.1339. 


lattice data set 

/(O) 

Ao/a^ 

A+/a^ 

0 = 0 

1.107(105) 

5.2(2.1) 

8.1(2.4) 

0 = 0 + kin. A, B 

1.089 (47) 

4.7(8) 

7.7(9) 

0 = 0 + kin. C 

0.9934(23) 

3.3(6) 

6 . 0 ( 8 ) 

0 = 0 + kin. A, B, C 

0.9966(25) 

3.0(5) 

5.5(7) 


Table 3; The same as in Table 2 but excluding in the data sets the value of /o(<?^) at 

H Simax' 


and the slopes Aq and A+, obtaining a precision competitive with that reported in 
Table 2, where the very precise value /o(9max) i® included. 

We have also performed the analysis of our lattice simulations at kg = 0.1339 and k^ = 
0.1349. The results obtained are similar to the ones shown in Figs. 2-4 and in Tables 2- 
3. Our findings about the impact of the twisted BCs remain unchanged. The same is 
true also for the K ^ K and tt ^ vr transitions, where g^ax vanishes, f-{q^) = 0 and 
/(O) = 1 because of vector current conservation. Note that for degenerate transitions, 
like the ones needed for the nucleon magnetic and the neutron electric dipole form factors, 
the kinematics A and B precisely coincide, while the kinematics C does not add any new 
information. 

We mention that there are other cases of phenomenological interest where a lat¬ 
tice point at zero-momentum transfer is not accessible with periodic BCs. For instance 
the nucleon magnetic form factor Gm{q‘^) is related to the matrix elements of the spa¬ 
tial components of the electromagnetic current, which are proportional to the value of 
the momentum transfer. Therefore, the nucleon magnetic moment Gm{q^ = 0) is not 
directly accessible with periodic BCs and it should be determined by a “long” extrapo¬ 
lation to = 0 (see Ref. [10]). Another example is the neutron electric dipole moment 
as determined with the strategies described in Refs. [11,12], which require again a long 
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extrapolation of the CP-violating neutron form factor F^^q'^) to ( 7 ^ = 0 In such cases 
however the introduction of twisted BCs is not trivial, because the current involved, the 
electromagnetic one, is not flavor changing. One can speculate that, by introducing an 
additional flavor and suitable interpolating fields, the application of the twisted BCs to 
the additional flavor, at least for quenched simulations, might provide the form factor of 
interest at quite small values of the four-momentum transfer. The application of twisted 
BCs to electromagnetic transitions requires therefore a careful treatment, which is well 
beyond the scope of the present work. 

6 Conclusions 

We have investigated the application of twisted boundary conditions to quenched lattice 
QCD simulations of three-point correlation functions in order to access spatial compo¬ 
nents of hadronic momenta different from the integer multiples of 27r/L. The vector 
and scalar form factors relevant to the K ^ n semileptonic decay have been evaluated 
by twisting in all possible ways one of the quark lines in the three-point functions. We 
have found that the momentum shift produced by the twisted boundary conditions does 
not introduce any additional noise and easily allows to determine within a few percent 
statistical accuracy the form factors at quite small values of the four-momentum trans¬ 
fer, which are not accessible when periodic boundary conditions are considered. We are 
confident that these findings are independent of the use of the quenched approximation, 
so that they will hold as well also in case of partially quenched and full QCD simulations. 

We have studied the impact of twisted boundary conditions on the precision of the 
determination of the K ^ n form factors and their slopes at zero momentum transfer. 
We have found that: i) when the precise lattice point for fo{q'^) at q'^ = is not 
included in the analysis, the use of twisted BCs is crucial to allow a determination of 
the form factor at zero-momentum transfer at a few percent statistical level; ii) when 
the very accurate value /o(9max) i® included in the analysis, the impact of twisted BCs 
turns out to be marginal, while we remind it is expensive for the computational time. 
The latter result is due to the fact that; i) the precision of the lattice points obtained 
with twisted BCs is not comparable to the one that can be achieved for /o((7max); ii) 
the use of twisted boundary conditions does not lead to a sufficient improvement of the 
precision with respect to the nearest points obtained with periodic boundary conditions. 

We stress that there are cases of phenomenological interest where a lattice point at 
zero-momentum transfer is not accessible with periodic BCs, like e.g. the case of the 
weak magnetism, the weak electricity, the induced scalar and pseudoscalar form factors 
entering hyperon semileptonic decays. 

direct determination of the form factor at = 0 can be obtained with the approach 

proposed in Ref. [13]. 
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